A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic axisymmetric response of viscoplastic shells of revolution. The analytical formulation is based on Reissner's large deformation shell theory with the effects of transverse shear deformation, rotatory inertia and moments turning around the normal to the middle surface included. The material model is chosen to be isothermal viscoplasticity, and an associated flow rule is used with a von Mises effective stress. A mixed formulation is used with the fundamental unknowns consisting of six stress resultants, three generalized displacements and three velocity components. Spatial discretization is performed using finite elements, with discontinuous stress resultants across element interfaces. The temporal integration is performed by using an explicit central difference scheme (leap-frog method) with an implicit constitutive update. The sensitivity coefficients are evaluated using a direct differentiation approach. Numerical results are presented for a spherical cap subjected to step loading, and a circular plate subjected to impulsive loading. The sensitivity coefficients are generated by evaluating the derivatives of the response quantities with respect to thickness, mass density, Young's modulus, and two of the material parameters characterizing the viscoplastic response. Time histories of the response and sensitivity coefficients are presented, along with spatial distributions of these quantities at selected times.